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Abstract. The Ha equivalent width (EW) is the ratio of the Ha flux to the con- 
tinuum at 6565A. In normal star forming galaxies the Ha flux is dominated by repro- 
cessed photons from stars with masses greater than 10 Mq and the 6565 A continuum 
is predominantly due to 0.7 - 3.0 Mq red giant stars. In these galaxies the Ha EW is 
effectively the ratio of high ma ss to low mass stars and i s thus sensitive to the stellar 
initial mass function (IMF). In i Hoversten & G lazebrookl ( l2008h we used ~ 131,000 
galaxies from the Sloan Digital Sky Survey to show evidence for systematic variations 
in the IMF with galaxy luminosity. In this proceeding we use that sample, with the ad- 
dition of H6a measurements, to investigate other parameterizations of the IMF. We find 
evidence for IMF variations with surface brightness, and also show that, modulo uncer- 
tainties in spectral synthesis models, that 120 Mq stars are important in accounting for 
the observed Ho- EW distribution. 



1. Introduction 

The stellar initial mass function (IMF) has historically been treated as universal across 
environments and time given that the large errors on the IM F slope from individual 
cluster measurements are all consistent with a single slope (iKroupalboOl ). Recent 



studies which infer the slope of the IMF from the integrated l ight properties of galaxies 



have provided evidence for systematic variations in the IMF. IHoversten & Glazebrook 



(12008,. hereafter HG08) found evidence for a link between galaxy luminosity and IMF 
slope based on a statistical analysis of the g - r colors and Ha equiv alent widths (EW ) 
of ~ 131,000 galaxies from the Sloan Digital Sky Survey (SDSS; Yorket al. 200(f). 
They found that the IMF slope was steeper in low luminosity galaxies. .Meurer et al] 



(I2009L hereafter M09) found that m a sample of HI selected galaxies the extinction 



corrected Ha to far-ultraviolet (FUV) flux ratio is correlated with galaxy luminosity 
which they attribute to IMF variations in agreement with HG08. However, they found 
that the correlation was stronger with galaxy surface brightn ess which they argue is the 
more physically fundamental connection. iLee et al.l (12009 ) compared star formation 
rates (SFR) calculated from FUV and Ha flux for a sample of galaxies spanning five 
orders of magnitude in SFR. They found that with decreasing galaxy luminosity the Ha 
SFR systematically underestimated the FUV SFR. One explanation for their results is 
a systematically varying IMF. 
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This proceeding expands on HG08, utilizing the same data set and techniques 
except where denoted. For an in depth description the reader is directed to HG08, 
however a brief overview is provided here. The HG08 sar nple contains 130,602 star 
formi ng galaxies with H/3 S/N > 5 selected from SDSS DR4 dAdelman-McCarthy et al. 



12006'). However Ha measurements, as well as the ¥15 a , were not from the standard 
SDSS pipeline, but rather calculated by iTremonti et all (12004 ). Values are extinction 
corrected and ^-corrected to z = 0. 1 . A large grid of galaxy spectral synthesis models 
were calculated over a range of ages, star formation histories (SFH), metallicities, and 



IMF slopes r using the publicly available PEGASE models (iFioc & Rocca-Volmerange 
19971) . 



HG08 compared the (g-r)n i c olor to the Ha EW, a method pioneered by lKennicutt 



( 1983 ) and iKennicutt et aL ( 1994 f). In this parameter space IMF variations dominate 



variations in metallicity and SFH time constant and are largely orthogonal to variations 
in extinction and age as shown in Figure 3 of HG08. The Ho- EW is the ratio of flux 
in the Ha line, which in star forming galaxies is due to reprocessed ionizing photons 
from young stars with masses greater than lOM©, and the continuum at 6565A which is 
dominated by red giant stars with masses in the 0.7-3 M© range. This is eff"ectively the 
ratio of the number of massive stars to those around a solar mass and is thus sensitive 
to the IMF. 

In HG08 the IMF is parameterized as 

dn f -0.5 for 0.1<m/Mo<0.5 



dlogm \ -r for 0.5 < m/Mg < 120 



following IB aldry & Glazebrookl (120031) . For a Salpeter-like IMF, which is surprisingly 
consistent with modern values, F = 1.35, however the original Salpeter IMF was a 
single power-law with no break. 

In this proceeding we consider two additional IMF parameterizations. As dis- 
cussed in HG08, in this color vs. Ha EW parameter space the IMF models themselves 
can be degenerate. In Equation [T] the IMF slope, F, is allowed to vary, while the upper 
mass cutoff is held fixed. Another option, described by 

dn f -0.5 for 0.1 < m/Mg < 0.5 



dlogm ^ \ -1.35 for 0.5 < m/Mo < M,p 

is to keep the slope at a Salpeter value (F = 1.35), but allow the upper mass cutoff, M„p, 
to float from 20 to 150 Mq. For M^p < 120 M© the models can be explicitly calculated 
with the pEGASE models. Above 120 M© the values are extrapolated. 

The second new IMF parameterization is the "minimal" case of the integrated 
galaxial initial mass function (IGIMF) theory ( Weidner & Kroupa 2005). In the IGIMF 
theory the stellar IMF is universal. However the combined effects of sampling the IMF 
from many star forming regions, which themselves have a size distribution, leads to 
an IGIMF which can differ from the universal IMF of stars. The IGIMF is dependent 
on the SFR of the galaxy as the maximum stellar cluster mass is related to the galaxy 
SFR. The minimal IMF case was ca lculate d for a range of SFRs using the pEGASE 
input files of Pflamm-Altenburg & Kroupal (l2009l) . A comparison of the three IMF 
parameterizations can be found in Figure [T] 

The remainder of this proceeding discusses two extensions of the HG08 study. In 
Section 2, the addition of H(5 absorption strength to the fitting procedure is discussed. 
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Figure 1. Model {g - r)o.i vs. Ha EW tracks for several IMF parameterizations. 
Each contour represents the area of parameter space covered by the 18,480 models 
(described in the text) for the given IMF. In the left panel the slope of the IMF is 
varied as in Equation [T] with an upper mass cutoff of 120 M©, in the middle panel 
r = 1.35 while the upper mass cutoff varies following Equation|2l and in the right 
panel the minimal IGIMF scenario is used. The hatched area covers the area of the 
parameter space for the F = 1 .35, M„p = 120 M© models and is identical in all three 
panels. 



Section 3 investigates IMF variations with galaxy surface brightness, followed by a 
discussion of the results in Section 4. 



2. H(5a Absorption 

Comparing the color to the Hor EW is a two parameter fit and thus has inherent limi- 
tations. Changing the age, metallicity, dust extinction, SFH, and IMF can have similar 
effects. Fortunately, in the {g - r)o.i vs. Ha EW plane used by HG08 the IMF slope 
largely distinguishes itself from the other effects. However in Figure 4 of HG08 it can 
be seen that for galaxies with red colors and low Ha EWs the different IMF model 
tracks become degenerate again. Of particular importance is the influence of bursts and 
gasps of star formation which can mimic many other effects. This motivates the use of 
a third parameter. 

Section 7 of HG08 explores the use of the H5 absorption feature to place further 
constraints on their results from {g - r)o.i vs. Ha EW modeling. H(5 absorption is due 
to stellar photospheric absorption lines. Balmer absorption lines are most pronounced 
in A stars, and weaken for both hotter and cooler stars following the Saha equation. 
For this reason the HJ absorption can be used as a proxy for the fraction of light being 
supplied by A stars, and less significantly by B and F stars, in a stellar population. In a 
simple stellar population H(5 absorption will peak at around 1 Gyr. At this time O and 
B stars have already burned out, but the A stars have yet to leave the main sequence. 
As such, H(5 absorption can be used to probe the age of a stellar population as well as 
uncover evidence of star formation bursts around 1 Gyr in the past. 
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IWorthev & Ottaviani (119971) set forth two different methods of measuring the 



absorption. As in HG08, we use the US a method which has a wider central bandpass to 
match the line profiles of A stars. This is the preferred index for galaxies because it is 
less noisy in low-S/N spectra and velocity dispersion widens absorptions features. The 
HSa index ranges from 13 for A4 dwarf stars to -9 in M-type giants. 

HG08 found that the distribution of US a was different for their highest and lowest 
luminosity galaxy bins. For the most luminous galaxies which are in good agreement 
with a universal Salpeter IMF the H5a values are centered at 3.4 and suggestive of 
a range of H(5a- On the other hand the distribution for the faintest luminosity bin is 
consistent with all galaxies having H5a=6.1. 

When trea ted separat ely, as in HG08, H5a provides only weak constraints on the 
IMF. However, iHoversteiil dlOOS) gave preliminary results demonstrating the power of 
including H^a explicitly in the (g - r)o.i vs. Ha EW analysis. In the three parameter 
analysis. Equation 8 of HG08 is replaced by 

;,?(r.z.,^) . - ^(r.z.>.^)j%|H.,. - Mr.z.,.»)j' ^ 0) 

in the "pseudo-;^f^" minimization, for a model with IMF F, metallicity Z, age t, and SFH 
(A. For each galaxy /, the inputs are the observed (g - r)o.i color c,. Ha EW w,-, and HSa 
index 5,- and measurement errors o"c;, cr„.. and cr^.. These are compared to the model 
values c{T, Z, t, ip), w(F, Z, t, i/f), and 6(T, Z, t, if/). From these the statistical estimator 

Z, t, ij/) value is calculated by brute force for each galaxy, and then mai^ginalized 
over metallicity, age, and SFH following HG08. 

The results of the three-parameter analysis are shown in Figure |2] The top row 
of Figure |2] shows the best-fitting IMF models as a function of the ro.i luminosity for 
each of the three IMF parameterizations. Each point represents 500 galaxies binned by 
luminosity, and the lines show the upper and lower 95% confidence limits on F. The 

bottom row shows the goodness of fit parameter, which is discussed in HG08 and 
should not be used for any conventional, textbook interpretation. 

Figure [2] shows that the inclusion of H(5a has a minimal effect on the derived F 
values at low luminosities. However, on the bright end inclusion of H(5a has a strong 
influence on the results. In the two-parameter fit of HG08, the IMF slope F steepens at 
the highest luminosities. In the three-parameter fit, F flattens out to nearly the Salpeter 
slope which much more closely matches what is expected in bright galaxies. The bot- 
tom panel of Figure[2]shows the quahty of fit parameter The fit to a single, universal 
IMF is much better for the high luminosity galaxies, and then deteriorates for the low 
luminosity galaxies which favor steeper IMF slopes. 

Figure |2] also shows evidence for an increase in Mup with increasing galaxy lumi- 
nosity when using the IMF parameterization of Figure|2] For the IGIMF model the SFR 
increases with galaxy luminosity as expected, however most galaxies are best fit with a 
SFR of 1000 M© which is a factor of 2 to 3 higher than what is generally accepted for 
these types of galaxies. 



3. Surface Brightness 



HG08 only investigated IMF variations correlated to galaxy luminosity. However, M09 
has argued that the surface brightness, via the surface mass density, is more fundamen- 
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Figure 2. MC simulation three parameter results for the full sample binned by 
M^- o.i. As in Figure [T] the left column gives results for the IMF of Equation [T] at 
center are the results assuming Equation|2l and at right are the results for the minimal 
IGIMF. Each diamond represents 500 galaxies plotted at the mean Mr^o.i value of the 
bin. The solid lines represent the upper and lower 95% confidence region measured 
for each bin. The three parameter results, including the HS absorption data, are 
plotted in black. In the top row the best fitting IMF model values as a function of 

r-band luminosity, while in the bottom row the;^'^ values for each luminosity bin are 
plotted. 



tally related to IMF variations than luminosity. Using the sample and techniques of 
HG08 it is straightforward to repeat the analysis replacing luminosity with the surface 
brightness. This provides an independent check on the surface brightness results in 
M09 using a completely different technique and galaxy sample. 

There are a number of different ways to measure surface brightness in SDSS galax- 
ies. Here we use the following definition 

j"r,50 ^r-Ar-kr + 2.5 log 2 + 2.5 log {nrj^) (4) 

where /i,-5o is the ro.i surface brightness within the Petrosian half light radius. The 
surface brightness is a function of the SDSS r-band Petrosian magnitude r, the r-band 



Milky Way extinction measured by Schleg el et al.l ( 1998), the r-band iir-correction to 
the z = 0.1 calculated from the code of Blant on et all (I2003D . and the r-band Petrosian 
half-light radius r^Q measured by the SDSS. 

The differences in the Ha EWs and colors for the highest and lowest surface bright- 
ness galaxies can be seen in Figure |5] The two panels show the 10,000 lowest (left), 
and 10,000 highest (right) fi^^so galaxies in the sample. The locus of the high surface 
brightness galaxies tracks the Salpeter IMF model. However, for fixed (g - r)o.i the 
locus of the low surface brightness galaxies is clearly lower in Ha EW, suggesting a 
steeper IMF slope. 
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Figure 3. Distribution in (g - r)o.i - log(HQ' EW) space for the 10,000 lowest (left 
panel) and 10,000 highest (right panel) fir.so galaxies in the sample. Contours are 
logarithmic. Outside the last contour individual galaxies are plotted. Solid lines 
are model tracks with exponentially decreasing SFHs with t = 1.1 Gyr and solar 
metallicity. The age increases along the tracks from 100 Myr in the upper left to 13 
Gyr at the lower right. The upper line has T = 1 .00, the middle line is similar to 
Salpeter's IMF with F - 1.35, and the lower line has F - 2.00. The cross in the 
lower left of each panel indicates the median error bars of the sample. 



In agreement with the empirical trends of Figure[3l Figure[3]shows the IMF results 
as a function of //^.so for the three-parameter analysis of Section 2. The relationship 
between iir,5Q and F is nearly linear before flattening out at high surface brightness. A 
similar relation is found between surface brightness and M„p. As in Figure|2j the IGIMF 
analysis yields SFR values which are implausibly high, although the qualitative trend 
of increasing SFR with increasing surface brightness is expected. In a departure from 
the luminosity based results in Figure |2] the quality of the fit is much more consistent 
across the surface brightness bins, particularly for the three-parameter fit. 

The SDSS spectroscopic fibers have a fixed size of 3 arcsec. Given our definition 
of the surface brightness, bright galaxies with pronounced bulges will have smaller 
Petrosian half-light radii and will have brighter /ir,50 values. This is why it is preferable 
to tailor the aperture to include as much light from the galaxy as possible, as in M09. 
The disadvantage is that it is a more labor intensive process; the M09 sample contains 
103 galaxies, while ours is over 1,000 time larger which enables more detailed statistical 
analyses. Given the nature of this sample disentangling the effects of surface brightness 
and aperture fraction is a challenge. 
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Figure 4. MC simulation three parameter results for the full sample binned by 
//r,50- The rest of the description is identical to Figure|2l 

In M09 it is argued that surface brightness are more fundamental than the lumi- 
nosity in IMF variations. Given our large sample size we can plot the best fitting F as a 
function of r-band luminosity and surface brightness simultaneously. This reveals that 
the variations in F are not orthogonal to either axis, showing that the IMF variations are 
a function of both parameters simultaneously. 



4. Discussion 

The analysis presented here shows that in large samples of galaxies there is a clear 
trend with both galaxy luminosity and surface brightness which cannot be explained 
using current spectral synthesis models with a universal IMF. However, there are two 
primary weaknesses of the results described above. 

The most significant limitation of the type of analysis presented here is the treat- 
ment of the SFH. The models described here assume a smoothly varying SFH. How- 
ever, as discussed in HG08, SFH discontinuities can have a large impact on the results. 
For instance, during a burst of star formation the Ha EW will increase and galaxy 
colors will become bluer which can change the best fitting IMF. Nonetheless, HG08 
demonstrates that while this is an intractable problem for individual galaxies, for statis- 
tical samples of galaxies it can be shown that reproducing the observed distribution of 
galaxies requires an implausible coordination of burst times. A new analysis of the SFH 
incorporating the ¥16a absorption and the new IMF parameterizations is being made for 
an upcoming paper. 

The second problem is that spectral synthesis models have inherent limitations 
and are still being actively developed. The systematic effects of the choice of spectral 
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synthesis models is not known for this type of analysis, and the results are subject to 
change based on improvements in models of stellar evolution. To quantify this effect, 
multiple models are being used to determine the scale of variations between models. 
The fundamental limitation that spectral synthesis models are only as good as the stellar 
models that go into them still remains. 

However, if variation of the IMF shown here is taken at face value, the results here 
provide insights into the IMF. Principle among these is that 120 M© stars are important 
for creating the observed Ha EWs. Figures [2] and [3] show that the most luminous and 
highest surface brightness galaxies favor the presence of 120 M© stars when assuming 
an IMF of the form of Equation [2l At the highest luminosities and surface brightnesses 
the fit is improved by allowing even more massive stars to form. 

Secondly, the qualitative behavior of the IGIMF parameterization is what one 
would expect. With increasing luminosity and surface brightness the galaxies favor 
IGIMFs with higher SFRs. However the calibration of the IGIMF seems to be off 
as most of the galaxies are best fit by an IGIMF model with a SFR of 1000 M© 
yr"^ which is largely unphysical. The "minimal" IGIMF model is truncated at 120 
M© and has more massive stars tha n the standard IGIMF model with 120 M© cutoff 
(iPflamm-Al tenburg & Kro"i^l2009h . ] lowever even in the minimal case at the highest 
SFR there are significantly fewer 120 M© stars than in an IMF from Equation [T] with 
r = 1.35 or from Equation |2] with M„p = 120©. While it needs to be explicitly mod- 
eled, it seems that increasing the maximum stellar mass in the IGIMF parameterization 
beyond M„p = 120© could solve this problem. 

The results here are also in agreement with M09, in that the IMF appears to vary 
as a function of galaxy surface brightness. Lastly, the addition of the h6A data in the 
analysis improves the constraints on the IMF slope for the most luminous galaxies. 

A more in depth analysis of these preliminary results is currently being performed. 
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